Integrability of zero-dimensional replica field theories at (3 = 1 
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Building on insights from the theory of integrable lattices, the integrability is claimed for nonlinear 
replica a models derived in the context of real symmetric random matrices. Specifically, the fermionic 
and the bosonic replica partition functions are proven to form a single (supersymmetric) Pfaff-KP 
hierarchy whose replica limit is shown to reproduce the celebrated nonperturbative formula for 
the density-density correlation function in the infinite-dimensional Gaussian Orthogonal Ensemble. 
Implications of the formalism outlined are briefiy discussed. 
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Introduction. — Perhaps the most comprehensive and 
unifying quantitative description of disordered and quan- 
tum chaotic systems can be achieved by using field- 
theoretic methods. Since the seventies of the last cen- 
tury, a number of such field-theoretic frameworks have 
been devised. Having different mathematical status and 
physical areas of applicability, they all rest on the concept 
of functional field integrals and are known as nonlinear 
a modelai"— . 

In the context of disordered systems, three major for- 
mulations of nonlinear a models exist: the replicai^— , the 
supersymmetry^, and the Keldysh^ a models. Of these, 
replicas is the earliest^i^ii and the most controversial in- 
vention that continues to challenge^"— both mathemat- 
ical and theoretical physicists. Decades after debuting 
in condensed matter physics, an operationally transpar- 
ent and controllable treatment of replica field theories 
is barely available-i even in the mathematically simplest 
yet physics motivated setup provided by the random ma- 
trix theorjiS, (RMT). Indeed, it took nearly twenty years 
to extract some nonperturbative results out of replicas, in 
both the RMT limiti^ and beyond i\^. More precisely, 
adopting the ideas of replica symmetry breaking origi- 
nally devised in the theory of spin glasses^, the authors 
of Ref.^'^ managed to show that the fermionic variation^ 
of replicas is capable of producing partially nonperturba- 
tive RMT results for the microscopic two-point density- 
density spectral correlation function 
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Equation ([l} describes two-point spectral correlations in 
'infinite-dimensional' real symmetric (/? = 1), complex 
Hermitean (/3 = 2), or quaternion- real self-dual (/3 — 4) 
random matrices whose spectrum was unfolded so as to 
make the mean level spacing unity, A = 1. On the formal 
level, R2 (io) corresponds to the thermodynamic limit 



Ri^\uj) = lim 

N^oc 



Rf\EuE2;N) 
Rf\Er,N)Rf\E2;N) 



(2) 



where Ei^2 = E± LuR^f\E; N)/2, and a set of finite-iV 
correlation functions is defined as follows: 
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Here, E stands for E = {Ei,--- ,Ep), whilst angular 
brackets (• • • ) denote averaging over the Gaussian or- 
thogonal (/3 — 1), unitary (/3 — 2), or symplectic (/? = 4) 
ensemble of random matrices specified by the normal 
probability density oc exp(— Tr D-C ) on the space MAr(F^) 
of Hermitean N x N matrices with real (Fi = K) , com- 
plex (F2 = C), or quaternion- real (F4 = H) entries. 

The validity of Eq. ([1]) , summarising achievements of 
the early fermionic replica research, is restricted^ to 
the energy scales in excess of the mean level spacing, 
w ^ A = 1. Notably, several attempts^'^ to reproduce 
the same result from bosonic replicas did not bear fruit. 
Yet, within the supersymmetry technique^, the same task 
of computing the density-density correlation function in 
the RMT setting was accomplished^ti^ in less than a year 
and with greater rigour. For the Dyson (3 — 1 sym- 
metry class (which will be a focus of this Letter), the 
truly nonperturbative density-density correlation func- 
tion readsi^iii (w > 0): 
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Sadly, heuristic approaches to nonlinear replica a models 
have so far failed to produce the result of such generality. 
One of the reasons behind this misfortune of replicas 
is that they are much more quality demanding in mak- 
ing various approximations and, consequently, are more 
involved operationally-wise. Such a vulnerability can be 
attributed to the continuous geometry underlying replica 
field theories. The latter derives from the very idea of 
constructing replica generating functionals which asso- 
ciates a physical observable of interest with the n — >■ 
limit of a collective matrix field Q arising as a result of 
imaginary n-fold replication of original physical system. 



Since the dimension of Q is an integer, dini(Q) oc n, tak- 
ing the rephca hmit becomes quite an adventure^i^iiiiiH. 
For instance, the p-point Green function 
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can formaUy be recovered from the rephca hmit 
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where Zn {E]N) is either fermionic (n > 0) or bosonic 
(n < 0) representation of the rephca partition function 
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The two representations can be derived from Eq. ([7]) 
by using standard held-theoretic methods^i^ii^. Alter- 
natively, one may appeal to the RMT duality relationsi^ 
which yield, e.g., the fermionic replica partition function 
{n > 0) in the form 



zI^^+\E;N) oc e(2//3)T.-[S8l„]M[^Q]g-(/3/2)TrQ = 

X det^Qe"2^Tr[Q(£«l„)]^ (g) 



Here, an integral runs over the matrix field Q G 
^np(]F4//3) of dimension dim(Q) = np. The bosonic 
replica partition function (n < 0) whose derivation is 
more intricate and delicate^^ also admits a dual matrix 
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akin to Eq. ^. 



The necessity of performing a limiting procedure in 
the matrix dimension [Eqs. ^ and dH)] highlights an 
unusual continuous geometry of replica field theories; it 
ruins a classic notion of matrix, raising both conceptual 
problems^- regarding mathematical foundations of non- 
linear replica a models and operational problems of deal- 
ing with weird objects where the plain intuition often 
refuses to work. 

To avoid running into an underdeveloped concept of 
spaces of rational dimensior^^^, one first performs replica 
calculations for an integer number of replicas in a hope to 
implement the replica limit n — >■ at a later stage after 
seeking an analytic continuation away from n integers. 
For such a continuation to be grounded, one should start 
with an exact integer-n result. The latter had been un- 
available in the early replica studies ^'^ that heavily relied 
on existence of a large spectral parameter uj/A ^ 1 re- 
quired to justify approximate saddle-point calculation of 
the replica partition function. The above inequality re- 
stricted domain of applicability of Eq. ([T]), while an ap- 
proximate character of integer-n calculations led to its 
questioned^ mathematical status. 

The status of nonlinear replica a models was dramat- 
ically uplifted in RefsJ2rJl where both fermionic and 



bosonic versions of zero-dimensional replica field theo- 
ries were shown to be exactly solvable for /? = 2 sym- 
metry class. The theory of Painleve transcendents^^, 
the Toda lattice-^- and associated r functions^S. lied 
at the heart of the exact approach to replicas in the 
elaboration3^^— . Later, a complementary - supersym- 
metric - formulation of replicas was introduced by Split- 
torff and Verbaarschot'^^. These authors have argued 
that the replica limit can efficiently be implemented 
on the level of supersymmetric Toda lattice equation 
whose positive and negative branches describe fermionic 
and bosonic partition functions, respectively. Super- 
symmetric replicas have greatly simplified calculations of 
/3 = 2 spectral correlation functions through a remark- 
able fermionic-bosonic factorisation22i^. 

In this Letter, we take one more step towards compre- 
hensive understanding of integrable structures of non- 
linear replica a models by presenting an integrable the- 
ory of a supersymmetric variation of zero-dimensional 
replica field theories for /3 = 1 symmetry class that so 
far denied a nonperturbative treatment. Focussing on 
the Gaussian orthogonal ensemble (GOE) and building 
on the theory of PfafF-KP r functiona^^, we shall show 
that supersymmetric replicas do produce a nonperturba- 
tive result Eq. ^ for the GOE two-point density-density 
correlation function. 

Density -density correlation Junction from supersym- 
metric replicas. — Specifying the GOEjv by the normal 
probability density oc exp(— A^TrJC^) on !MAr(Fi), we 
look for a large-A^ limit of the replica partition func- 
tion Zn {Ei,E2]N) whose energy variables £'1^2 are 
rescaled as described beneath Eq. ^. A sometwhat 
lengthy calculations! based on the methods described in 
Refj ^^'^^1'^^ shows that both fermionic and bosonic parti- 
tion functions admit the large-TV factorization 

z\^h^;N) = iv2"^(4e)Tl"l^ . z(±)(c.), (9) 
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where z\n (^) is the fermionic partition function, 
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Here, cj M> cj + iO is assumed to have an infinitesimally 
small positive imaginary part, A^^ = YYjyki'^j ^ ^fe) is 
the Vandermonde determinant, 
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In order to recover the two-point density-density cor- 
relation function [Eq. ^] in infinite-dimensional GOE, 
i?2(w) = (I/27r^) y\eg{uj) - 1/2, through the replica limit 
for the corresponding two-point Green function 
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one has to find a proper nonperturbative representation 
of either fermionic^^ or bosonic^^ partition functions, or 
opt for supersymmetric replicas?^. It is the latter route 
that will be explored below. 



To make a forthcoming presentation transparent, we 
first quote our main analytic result. Let z„(aj) be a su- 
persymmetric replica partition function, 
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For all n integers (both positive and negative), the follow- 
ing nonlinear recursive differential equation of the Pfaff- 
KP type holds: 
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Equation p6|) . whose derivation will be sketched below, 
is central to implementing the replica limit Eq. (J14p . In- 
deed, assuming that Eq. p^ stays valid for n real, we 
introduce a small-n expansion 



logz„(w) = nh{Lj) + r? h{io) + 0(n3) 



(17) 



and substitute it into Eq. (fT6|) to realize that /"(i^) 
should be set to zerc'^"' to ensure existence of the limit 
Eq. (J14p . Further, we make use of the boundary condi- 
tions for z[~}{u:) at cj — >■ 00 [see Eq. (ITT|) ] to figure out 
that /{(O) = -ITT. Finally, we combine Eqs. dUl) and (fT7|) 
to observe the relation g'{ijj) = —f2{uj) and subsequently 
derive from Eq. (IK 
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Hence, the Green function can solely be expressed in 
terms of replica partition function for one fermionic 
[z{ (w)] and one bosonic [£}~'^(a;)] flavor. Such a 'factor- 
ization property' was first observed in Refi^ in the con- 
text of supersymmetric replicas for j3 = 2 matrix models. 
To complete the evaluation of i?2(w), we use Eqs. P^ 
and (jlip to find out 

z^^\u:) ^ -^5'(c.), (19) 
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where S{uj) — sin(7raj)/(7ra;). Equations ([T5)) - (PU)). im- 
ply 
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^(s\u:) + S'{u:)p dtS{t)\ (21) 



or, equivalently, 

R2{uj) - i?2(oo) = ,5(0;) - (s'^iuj) + S'{uj) I dt S{t) 

(22) 
Setting R2{oo) to unity concludes our replica deriva- 
tion of the celebrated nonperturbative formula [Eq. ^] 
for the two-point density-density correlation function in 
infinite-dimensional GOE. 

Supersymmetric Pfaff~KP equation. — Now we sketch a 
derivation of the Pfaff-KP equation [Eq. P^ ]. In doing 
so, we follow a r function approach exposed in Refi^ (see 
also a tutorial paper—). We proceed in four steps. 

(i) — First, we define the fermionic and the bosonic t 
functions given by 

4t^(^;*)- An / dA.(l-AI)e-^^+2^(*^^^)|A^J^ 

(23) 
and ($He s > 0) 



^t.'<-')^(S;i^Il/ 



dXk 



V4^ 



sAfc-V(t;Afc) I A A I 

(24) 



respectively. Here, s = —iTruj/2, whilst V{t; A) = 
^- j^tjA-' is a deformation potential parameterized by 
infinitely many parameters t — {ti,t2, ■ ■ ■)■ To make con- 
nection with the supersymmetric replica partition func- 
tion z„(w), we also construct the supersymmetric r func- 
tion 



'^2imi(*;*)' ™ < 0; 



T2m(s;<) 



'2|mp 



1, 



TO = 0; 



(25) 



such that the following projection relation holds: 

Zn{uj) = f2n{s = ~iTTUj/2; 0). (26) 



(ii) — Second, we utilize the approach by Adler and van 
Moerbeke22, to claim that both fermionic and bosonic r 
functions satisfy an infinite set of hierarchically struc- 
tured nonlinear differential Pfaff-KP equations in the 
variable t. The first two members of the Pfaff-KP hi- 
erarchy, written in terms of Hirota derivatives-^, read: 

{Dt+3Dl-4D,Ds)rtKrtJ ^2^4^1,04^1, (27) 

and 
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Taken together with the identity^i^ Di T2m+2° T2m-2 = 
—Di T2m_2 T2m+2, the abovc implies that the supersym- 
metric r function T2m{s;t) is governed by a single pair 
of Pfaff-KP equations obtained from Eqs. (gT]) and (^5)) 
by dropping (±) everywhere. Deciphering Hirota deriva- 
tives, one comes down to 

{KP,Ulogr2„Mt) = i2lH!!i4l4%^±M, (29) 



and 
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where (KPi)t and (KP2)t are the first and the second 
Kadomtsev-Petviashvili operators given bj^ 
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(Hi) — Third, given the invariance of f,^' under the 
change of integration variables Afc t-^ fik+e{p1 — i-) fJ-l^ in 
Eqs. ([23l) and (|24|) . we observe-^ that the supersymmetric 
T function Eq. ([^5]) satisfies an infinite set of Virasoro 
constraints {q > —1) 
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dtq+1 , 
T2m(s;t)=0, (33) 



5t,at„_, ^('^+1) J: (^^) 



is the /? = 1 Virasoro operator—. Equation ([33]) assumes 
that djdtx = 9/9s; the operator 9/i9io should be inter- 
preted as djdta = 2m. 

(iv) — Fourth, we combine Eq. (|26|) with the three 
lowest Virasoro constraints [q — —1, 0, and +1) to 
project [(KPi)t-F(s/2n)(KP2)t] logf2„(s;t) onto the 
hyperplane t = 0. Lengthy but straightforward cal- 
culations yield^i the sought Pfaff-KP equation for the 
supersymmetric replica partition function Zn{uj). This 
completes our derivation of Eq. (ITC)) . 

Conclusions. — In this Letter, we have shown how the 
ideas of integrability can be utilized to formulate a non- 
perturbative theory of /3 = 1 zero-dimensional replicas 
in their supersymmetric elaboration. Although a partic- 
ular emphasis was placed onto the GOE random matri- 
ces, the formalism outlined is quite general and should 
equally apply to other random matrix models (including 
those appearing in various RMT formulations of quantum 
chromodynamics) , for both /3 = 1 and j3 = A symmetry 
classes which are dual to each other. Yet, we cannot 
but say that an understanding of formal structures lurk- 
ing behind zero-dimensional replica a models that was 
accumulated during the past decade takes us one step 
closer to formulating exact replica theories for more real- 
istic matrix models (e.g., random banded matrices) and 
physical systems. 
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